X-ray reflected spectra from accretion disk models. 
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: ABSTRACT 
O ■ 

■ We present new models for illuminated accretion disks, their structure and repro- 
^ ■ cessed emission. We consider the effects of incident X-rays on the surface of an accretion 

'— disk by solving simultaneously the equations of radiative transfer, energy balance and 

, ionization equilibrium over a large range of column densities. We assume plane-parallel 

\ geometry and azimuthal symmetry, such that each calculation corresponds to a ring at a 

QQ \ given distance from the central object. Our models include recent and complete atomic 

I data for K-shell processes of the iron and oxygen isonuclear sequences. We examine the 

' effect on the spectrum of fluorescent Ka line emission and absorption in the emitted 

■ spectrum. We also explore the dependence of the spectrum on the strength of the inci- 
. dent X-rays and other input parameters, and discuss the importance of Comptonization 

• • , on the emitted spectrum. 

. yr , 1. Introduction 

The X-ray spectra from active galactic nuclei (AGN) and X-ray binaries often show evidence 
of interaction between radiation emitted near the compact object and the nearby gas, which leads 
to signatures imprinted on the observed spectrum. The effects of this reflection include iron K 
line emission, in the range 6-8 keV, which is observed from nearly all accreting compact sources 



(jGottwald et al.lll995l ): and the flattening of the spectrum above 10 keV, usually called high energy 
bump or the Compton hump, since it originates due to the Compton scattering of photons by cold 
electrons in the gas. 



The profile of the iron K line carries important information about the physics of the material 
around the compact object. Since the emission may occur in a relatively small region close to 
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the center, the hne profile can be affected by relativistic effects. The be st-known example is the 
Seyfert 1 galaxy MCG-6-30-15 (jTanaka et al.lll995l : llwasawa et alJll999l ). for which the iron line 
appears to be broad and skewed well beyond the instrum ental resolutio n . Mo del li ne pro f iles fo r 
Schwarzschil d and Kerr metrics hav e been calculated by iFabian et al.l (|l989l ) and iLaoii (|l99ll ). 
respectively. iBrenneman &: Revnoldsl (|2006l ) have used this information to constrain the spin of the 
black hole harbored by this AGN. 

Observations of other accretin g systems such as low mass X-ray bi naries (LMXB) have also 
revealed the iron K line in emission. iBhattacharvva &: Strohmaverl (|2007l ) reported the detection of 
a broad iron line emi ssion in the XMM Ne wton spectrum of the LMXB system Serpens X-1. Using 
Suzaku observations, ICackett et al.l (j2008l ) confirmed this detection and found similar emission in 
two other LMXBs (4U 1820-30 and GX 349+2). Assuming relativistic broadening and applying 
the same models used in AGN observations, these authors have derived estimates for the inner 
radius of the accretion disk which can be used as an upper limit for th e radius of the neutr on star. 
Similar behavior appears in the millisecond pulsar SAX J1808. 4-3658 (jCackett et al.ll2009l ). Since 
the temperature of the surface layers of an accretion disk around a stellar-mass black hole or a 
neutron star is expected to be much higher than the temperature in the accretion disk of an AGN, 
the Compton down scattering of phot ons can also modify t h e emission profile and even contribute 
to the broadening of the iron K line (jRoss FabianI 120071 ). iReis et al.l (j2009l ) found evidence for 
such a situation in the Suzaku spectrum of the LMXB 4U 1705-44. 

An X-ray line propagating in a dense gas has a non-negligible probability of interacting with 
an electron, which leads to a down scattering of photons. The maximum energy shift per scattering 
of a photon with energy Ec due to the electron recoil (for backscattered photons at 180°), has 
magnitude AEmax = 2E'^/{meC^ + 2£'c)) where mgC^ is the electron rest-mass energy. Many 
scatterings will therefor e produce a discernibl e "Compton shou lder" between Ec and Ec — AE^ax 
(jPozdniakov et al.lll979l : llllarionov et al.lll979l : iMatt et al.lll99ll ). Although the Compton shoulder 
and the iron K/3 emission line are weaker than the iron Ka emission line, b oth have been detected 



in AGN observatio i is, especially those from the X-ray observatory Suzaku (IMarkqwitz et al. 



2007 



Reeves et al.ll2007l : lYaqoob et al.ll2007l ). and at least in one LMXB, GX 301-2 (IWatanabe et al, 
20031 ^. 



Initial studies of repro cessed radiation by cold matter were due to iLightman Sz Rvbickil (jl980l ) 
and ILightman et al.l (|l98ll ). who derived the Green's functions for th e scattering of photons b y 
cold electrons ari d disc ussed the implications for AGN observations (jLightman White! Il988l ) . 
George &: FabianI (jl99ll ) included line production in the ir Monte Carlo calculations, while pho- 
toionization equilibrium was included by lZycki et al.l (119941 ). Sever al authors expanded these studies 



1993 : lMatt et al 
1996; Matt et al. 



11996 



Poutanen et al.lll996l : iBlackman 



Krolik et a. 



static calculations hav e been carried out by 



Navakshin KallmanI mOl\ ): iBallantvne et al.1 (|200lh : 



19941 : 



I999I: 



assuming constant density along the vertical direction of the disk (iDone et al 
1993l:ICzernv fc Zvckll994 



1992 



Ross &: Fabian 



Rozanska Sz Czern' 



Dumont et al 



Magdziarz Sz Zdziarskilll995l : lRoss et al 
Ross Fabiarj 2005), while hydro 



Navakshin et al 



(2001 



((2003) and IRoss k Fabian 
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(120071). 



Nevertheless, one of the most significant hmitations inherent in current reflection models is 
the limited treatment of the physics governing the atomic processes affecting the excitation and 
emission from the ions in the atmosphere, especially the iron K-shell structure. In most cases, re- 
flection mo dels implement ana l ytic fi ts to the partial Hartree-Dirac-Slater photoionization cross 
sections of IVerner YakovlevI (119951 ). together with the line energies, fluorescence and Auger 
vields from iKaastra &: Mewd (19931) . which have been proven to be over simplified and incom- 



plete (jGorczvca k. McLaughlin 



20051 ). In some cases, the models lack atomic data for any of the 



neutral species (except for H and He). 

In this paper, we present new models for illuminated accretion disks a nd their structure im - 
plementing s tate-of-the-art atomic data for the isonuclear sequences of iron (iKallman et al.l l2004l ) , 
and oxygen (jGarcia et al.ll2005l ). In addition, the energy, spatial and angular resolution of our 
calculations are greater than previous works. The ionization balance calculations are performed in 
detail using the latest version of the photoionization code xstar. The radiation transfer equation is 
solved at each depth, energy and angle using the Feautrier formalism. We assume plane-parallel ge- 
ometry and azimuthal symmetry. For simplicity in this paper we assume constant density along the 
vertical direction of the atmosphere, and we leave the hydrostatic calculation for future publication. 

This paper is organized in the following way. In § [2] we describe the theory and numerical 
methods used to solve radiation transfer, ionization and energy equilibrium, and the atomic data. 
In §[3] we present a set of models for constant density atmospheres at different degrees of ionization, 
viewing and incidence angles, as well as the effect of the abundance on the reflected spectrum. 
Finally, the main conclusions are summarized in § [H 



2. Methodology 

In this section we describe the theory used in our calculations, in particular, the numerical 
methods to solve the radiation transfer equation, the ionization equilibrium, energy conservation 
and the new atomic data incorporated in our models. 



2.1. Radiative Transfer 

The standard form of the transfer equation for a one dimensional, plane-parallel atmosphere 
is given by 

^ar(^^ = viz, E) - x{z, E)Iiz, f,, E), (1) 

where /i is the cosine of the angle with respect to the normal, and ry(z, E) and xi^? E) are the total 
emissivity and opacity, respectively. It is convenient to write this equation in terms of the Thomson 
optical depth, dr = —axdz = —axnedz, where ax is the Thomson cross section (= 6.65 x 10~^^ 
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where we have defined 



cm ), and rig is the electron number density, such that 

^ir.E).^. (3) 

To solve the equati on of radiative transfer we use Lambda iteration in the Feautrier formalism, as 
described in § 6.3 of lMihalasI (|l978l ). By restricting /i to the half-range < < 1 to express the ra- 
diation field in its incoming and outgoing components, and using the symmetric and antisymmetric 
averages: 

u{T,fi,E) = ^ [/(r, E) + I{r, -ti,E)] (4) 

v{r, ^, E) = i[/(T, E) - I{t, -fi, E)], (5) 
one can rewrite the Equation ([2]) as a second order differential equation 

,WiEf^^ + ,MEf-^'-^ = u[,,E) - S[E) (6) 

where we have omitted the obvious dependence on r. Notice that the second term on the left hand 
side of the equation appears since we are using the energy-independent Thomson optical depth as 
length variable, instead of the total optical depth dr^E) = —x{z, E)dz. 

The second term in the right-hand side is the source function, which is defined at each depth 
and frequency as 

S{E) = ^ = ''knJ{E)+j{E) 

X{E) akn + Q-a 

where ^ 

J{E) = [ u{i^,E)dfi, (8) 
Jo 

is the mean intensity (first moment of the intensity), and j{E) is the thermal continuum plus line 
emissivity calculated at each depth. Here Okn and Oa are the scattering and absorption coefficients, 
respectively: the former is defined by the product of the density times the Klein-Nishina cross 
section, while the latter contains the continuum cross section due bound-bound, bound-free and 
free- free absorption. Both j{E) and are obtained from the ionization balance at each point using 
XSTAR. 

We also consider the redistribution of the photons due to Compton scattering using a Gaussian 
kernel convolved with J{E). The Gaussian is centered at 

S, = ^o(l + 40-eo), (9) 
where 6 = kT /m(.(? is the dimensionless temperature, Eq is the initial photon energy and eg = 

1 /2 

-Eo/"^eC^• The energy dispersion is given hy a = Eq \20 + |eg] 
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Therefore the total source function at each depth, as a function of the energy, is given by: 

S{E) 

where 



ME) + '^^^ 



X{E) '•^ ' x{E) 
J^{E) = j dE'J{E')P{E',E) 



(10) 
(11) 



is the Comptonized mean intensity, for which the kernel function (normalized to unity), can be 
written explicitly as 



-P(-E'c, Eg 



1 



crvr 



1/2 



exp 



-{Eg — Ec 



(12) 



As shown by IRoss FabianI (jl993l ). this treatment describes adequately the down sc attering of 
photo ns wi th energies less than ~ 200 keV. This approximation is also discussed by iRoss et al 



(j 19781 ) and iNavakshin et al.l (|200d ). There is an important limitation while using this treatment: 
if the bin width of a certain energy E^ is greater than 2 times the dispersion of the Gaussian, then 
there is a large probability for a photon to scatter into its own energy bin, producing a numerical 
pile up of photons at those energies. We use a logarithmically spaced energy grid, so the resolving 
power 7^ = E/ /S.E is a constant that only depends on the number of grid points. Therefore, the 
bin width is /S.E = E/TZ, and the pile up occurs at a critical energy Ep when /S.E = 2a{Ep); or 



Ep = TTleC 



47^2 



5kT 



1/2 



(13) 



From the last equation it is clear that the pile up appears only at low temperatures, specifically for 
temperatures less than Tp = nieC^ / ^kTZ^ . The lowest resolving power used in this paper isTZ = 350, 
and then Tp ~ 10^ °K. We will show later that we consider no cases for which the gas temperature 
becomes lower than this value. In fact, in the extreme limit of T = (and using the same resolving 
power), the highest pile up energy would be Ep ~ 1.6 keV, leaving the high energy part of the 
spectrum unaffected. 

To complete this solution, two boundary conditions are imposed; one at the top (r = 0), and 
one at the maximum depth (r = Tmax)- At the top, the incoming radiation field 1(0, — /i, E) is equal 
to a power law, with photon index and normalization as free (input) parameters in the calculation. 
Subtracting Equations @ and ^ 



uiO,n,E)-v{0,ii,E) = Ii, 



and it is easy to show that 
therefore at the surface, 



v{t,h,E) = u}{t,E)h 



du{T, n, E) 
97 



Uj{<d,E)lJL 



du{T, fi, E) 



u{0,fi,E) = -li, 



(14) 



(15) 



(16) 
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At the lower boundary we specify the outgoing I(r, +;U, E) radiation field to be equal to a blackbody 
with the expected temperature for the disk: 



du{T, ^, E) 



+ u{Tmax,fJ',E) = B(T, 



disk) 



(17) 



This condit ion appears since we assum e that there is an intrinsic disk radiation due to the viscosity 
of the gas (IShakura Sz Sunyaevlll973l ). The disk temperature at the lower boundary is then given 
by 



_j 3GMM \^ fRo 
''''' ~ \ 87raR^ [ [r 



1/2 



1/4 



(18) 



where G is the Newtonian gravitation constant, a is Stefan's constant, M is the mass of the central 
object, M the mass accretion rate, R is the distance frorn the ce nter, and Rq is the smallest radius 
at which the disk dissipates energy. IShakura &: Sunyaev predicted that for a non-rotating 

black hole this occurs at the innermost stable circular orbit (ISCO), or Rq = SRg = 6G M/c^, where 



■R.s is the Schwarzschild radius. However, recent magneto hydrodynamics calculations (jNoble et al 



20 id ) , have shown that the electromagnetic stress responsible for the energy generation rises steadily 
inward in the region inside the ISCO, falling sharply to zero just before the event horizon. Therefore, 
we will assume Rq = Rs, allowing energy generation in the inner part of the disk. It is clear that the 
properties of the accretion disk are introduced in the calculation by means of the lower boundary 
condition. Further, since we assume constant density for the gas along the vertical direction, it is 
customary to parametrize each model by the ratio of the net flux incident at the su rface over the 
density, using the common definition of the ionization parameter (jTarter et al.lll969l ): 



rip. 



(19) 



where i^x is the flux of the illuminating radiation in the 1-1000 Ry energy band. Finally, Equa- 
tions ([6]). (fT0|) . (fT6]) and (fT7|) are converted to a set of difference equations by discretization of depths, 
energies and angles. The solution of the system is found by forward elimination and back substi- 
tution. A full transfer solution requires the Feautrier solution to be iteratively repeated, in order 



to self-consistently treat the scattering process. This procedure requires 
iterations) for convergence. 



iterations (lambda 



2.2. Structure of the gas 



Given the solut ion for the radiation field a t each point in the atmosphere, we use the photoion- 
ization code XSTAR (jKallman &: Bautistall200ll ) to determine the state of the gas for a given value of 
the number density. The state of the gas is defined by its temperature and the level populations of 
the ions. The relative abundances of the ions of a given element and the level populations are found 
by solving the ionization equilibrium equations under the assumption of local balance, subject to 
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the constrain of particle number conservation for each element. Schematically, for each level 

Rate in = Rate out (20) 

The processes include spontaneous decay, photoionization, charge transfer, electron collisions, ra- 
diative and dielectronic recombination. Similarly, the temperature of the gas is found by solving 
the equation of thermal equilibrium, which may be written schematically as 

Heating = Cooling (21) 

The heating term includes photoionization heating (including the Auger effect), Compton heating, 
charge transfer and collisional de-excitation. The cooling term includes radiative and dielectronic 
recombination, bremsstrahlung, collisional ionization, collisional excitation of bound levels and 
charge transfer. All processes include their respective inverses so that the populations approach to 
local thermodynamic equilibrium (LTE) values under the appropriate conditions (i.e., high density 
or Planckian radiation field). 



2.3. Radiative equilibrium 

XSTAR calculates level populations, temperature, opacity and emissivity of the gas assuming 
that all the physical processes mentioned in the previous section are in steady state. Radiative 
equilibrium is achieved by calculating the integral over the net emitted and absorbed energies in 
the radiation field {Ec and E^, respectively); and varying the gas temperature until the integrals 
satisfy the criterion 

fi— 1£ < 10-^ (22) 
Eh + Ec 

The radiative equilibrium con dition must be ensured also while solving the transfer equation. Ac- 



cording to equation (2-83b) in iMihalasI (|l978l ): 



/■oo 

/ x{E)[JiE)-SmdE = 0, (23) 
Jo 

i.e., the total energy absorbed by a volume of material must be equal to the total energy emitted 
(note that, as in ^ 12.11 we have omitted the explicit dependence on r in the equations). By using 
the definitions of the flux and radiation pressure (second and third moments of the radiation field) : 

H{E) = [ v{fi,E)ndfi (24) 
Jo 

K{E) = I u{fi,E)fi'^dn, (25) 



and 







one can rewrite Equation (llSp as 



HiE)=.iE)^. (26) 
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Taking the derivative with respect to r and multiplying by u}{E) 



and comparing with the transfer Equation ^ integrated over 



i^{E)^^^ = J{E)-S{E). (28) 



From here it is clear that the radiative equilibrium condition (j23p is equivalent to 

— j^ H{E)dE = 0, (29) 

i.e., the net flux must be conserved for any optical depth. Further, inserting the source function (jlOp 
in the right-hand side of (I28p : 



Q poo roo 

— H{E)dE = — {akn[J{E)- UE)] + aaJ{E)-j{E)}dE. (30) 
OT Jq ar Jo 



If the gas temperature is high enough (T > 10^ °K), the Compton scattering is the dominant 
process, while both the opacity and the emissivity of the gas are negligible and only the first term 
in the right-hand side of Equation ()30p is important. For lower temperatures the inverse occurs: 
J{E) ~ Jc{E), which cancels out the first term. Thus, using (j29p 



1 

— [aaJ{E)-j{E)]dE = 0. (31) 

The conservation of the flux through the atmosphere is formally equivalent to the thermal equi- 
librium condition (I2ip . but its accuracy will depend on the error associated with the evaluation 
of the left-hand side of Equation (|3ip . Although we require a small error in the calculation of 
the emissivities j{E) and opacities aa by imposing the condition (j22p . the numerical conservation 
of the flux is a difficult task due to the fact that small errors accumulate over the large column 
densities covered in the calculation. It can be seen from Equation (j3ip that the error required in 
the calculation of emissivities and opacities must be of the order of a^^. Assuming that typical 
values for the density are rig ~ cm ^, then 10^ < arp^ < 10^^, much greater than the 

value required in ()22p . Therefore, we also apply a renormalization of the emissivities calculated by 
XSTAR before the solution of the transfer equation, such that: 

aaJ(E)dE 

r'"iE)=j{E) ^' I , (32) 
Jo jiE)dE 

at each depth. By implementing this correction in our models, we have found that the net flux (and 
therefore the energy) is conserved better than 1% for intermediate to large values of the ionization 
parameter (log ^ > 3) , but the maximum error can be ^ 7% for low- ionization calculations (log ^ < 
3). These errors are estimated by comparing the total fluxes incident on the gas (i.e., + Fdisk), 
with the total emergent flux. 
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2.4. Atomic data 



The XSTAR atomi c datab ase collects rec ent dat a, from many sources including CHIANTI 
dLandi k PhilliDslbood l. AD AS dSummerslbood'l. NIST dRalchenko et alJbood ). TOPbas e dCunto et al 



19931) and the IRON project (jHummer et al.lll993l ). The database is described in detail bv lBautista &: Kallman 



(|200ll ). Additionally, the atomic data associated with the K-shell of the Fe ions incorporated in 
the current version of XSTAR has been recently calculated and represents the most accurate and 
complete set available to the prese nt. Energy lev e ls and transition probabilities for first row ions 
Fe XVIII- Fe xxiii were reported b y iPalmeri et alJ (l2003al ) ; for second row ions Fe x-Fe xyii wer e 
reported by lMendoza et alJ (120041 ): and for the third row ions Fe ii-Fe ix by lPalmeri et al.1 (|2003bl ). 
The impact of th e damping by s pecta tor Auger resonances on the photoionization cross sections 
was discussed by iPalmeri et al. d2002l'l. P hot oi o niz at i on and electron impact cross sections were 
presented for second row ions bv lBautista et al.l (|2004l ) . Energ y levels, transition pro babilities and 
photoionization cross sections for Fe xxiv were calculated by iBautista et al.l (|2003l l . A compila- 
tion of t hese results an d a ca reful study of their impact on the photoionization models can be 
found in iKallman et al.l (120041 ) . Moreover, xstar also includes the atomic da.ta re levant to the 
photoabsorption near the K edge of all oxygen ions calculated by Garcia et al. (j2005l ). 



The approach used by these authors for the computation of the atomic parameters of iron 
and oxygen was based on the implementation and compa rison of result s obtai ned from dif ferent 
atomic codes of public domain, namely a utostructure (Badnell 1997 ) , HFR ( CowanlE98ll ). and 
the Breit-Pauli i?-Matrix package (bprm. ISeatonlll987l : iBerrington et al.lll987l ). 



2.5. Iteration procedure 

Starting at the top of the disk, the vertical structure of the gas is found by solving ionization 
and thermal balance at each spatial zone using XSTAR, as described in § 12.21 Since the radiation 
field is unknown the first time this is done, it is set to be equal to the incident power law at each 
depth. Once the last zone is reached, the temperature and density profiles are known, as well as the 
emissivities and opacities for each depth and energy. With this information the radiative transfer 
Equation ([6]) is solved as described in § 12.11 until the solution converges, which requires ~ T^^ax 
Lambda iterations. This provides a new and more accurate radiation field, that can be used to 
recalculate the structure of the gas. Because the emissivities and opacities are updated when the 
structure of the gas is recalculated, the radiative transfer calculations must be also repeated. We 
then continue with this process until all quantities stop changing within a small fraction. Typically, 
this requires ~ 20 gas structure calculations, times t^^x Lambda iterations. 

For our purposes, we desire to cover the largest optical depth possible for realistic computa- 
tional times and resources. Since Lambda iteration requires T^ax loops to converge, we performed 
different nu merical experimentatio n changing the value of Tmax- Comparisons with the semi-analytic 



solution by IChandrasekharl (jl960l ) showed that the error behaves, in general. 
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spending to a leakage of energy through a scattermg dominated slab. The comparison also showed 
that diminishing improvement is achieved after Tmax ~ 10 — 20, taking into account the required 
number of iterations. Therefore we limit all our calculations to Tmax = 10, with the understanding 
that our resulting global energy budget may be in error by as much as 10%. Finally, our calculations 
consider high resolution spectra with an energy grid of at least 5 x 10^ points (7^ = E/AE ~ 350), 
200 spatial zones, and 50 angles to account for anisotropy of the radiation field. Typically, each one 
of these models can be calculated in few hours in a last generation PC, which means that a grid 
of 20-30 models can be produced in one or two weeks, depending on the computational resources 
available, and the specific parameters to represent the models (e.g., lower ionization models tend 
to run slower than high ionization cases). This is relevant since XSPEC model tables can be easily 
generated and used in the interpretation of astronomical observations. Higher resolution can be 
achieved, although computational time increases (at least) with the square of the number of energy 
grid points, quickly imposing limits on the sizes of the grids. 



3. Results 

In this section we show the results obtained with our model for constant density atmospheres. 
We present a total of 20 reflection models, calculated for various representative conditions. These 
models and input parameters used to produce them are summarized in Table [H The first column of 
the Table contains an identification number assigned to each model and it will be used for reference 
through the rest of this paper. The next columns contain the input parameters that have been 
varied for each model, namely: the flux of the ionizing radiation (in the 1-1000 Ry energy 
range), the resulting logarithm of the ionization parameter ^ (see Equation [T9|), the spectral energy 
resolution or resolving power TZ = E/AE (determined by the number of energy grid-points used), 
the cosine of the incidence angle of the illuminating radiation with respect to the normal of the 
disk fiQ, and the abundance of iron normalized to its solar value Ape. Other input parameters not 
listed in the Table since they are constant and common to all the models are: density rig = 10^^ 
cm~^, photon index of the incident radiation F = 2, mass of the central object M = 10^ M©, 
distance from the central object R = 7Rs, and the mass accretion rate M = 1.6 x MEdd, where 
Rs = 2GM/(P' is the Schwarzschild radius and M^dd is the accretion rate at the Eddington limit. 
We have chosen these parameters such that the numbers are similar to those typically used by 
previous authors. However, note that these parameters only affect the present models by changing 
the effective temperature of the disk (equivalent to changing the intrinsic disk flux), to be used 
as the inner boundary condition (jl7p in the radiation transfer problem. Specifically, by using 
Equation (llSh this set of parameters correspond to an intrinsic disk flux of Fdisk = 3.6 x 10^^ 
erg/cm^/s, or an effective temperature of Tdisk = 2.8 x 10^ °K, which represents a cold disk when 
compared to the temperatures typically produced by the illuminating fluxes used in this paper. 
This is convenient since it allows us to analyze the reprocessed spectrum as a direct consequence of 
the incident power law without significant modifications due to the black body of the disk. 
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Table 1. List of reflection models with their respective input parameters 



Model 




log 


n 


Mo 


Ape 


1 


5 X 10" 


0.8 


350 


0, 


.71 


1.0 


2 


1 X 10^5 


1.1 


350 


0, 


.71 


1.0 


3 


2 X lO^-"^ 


1.5 


350 


0, 


.71 


1.0 


4 


5 X 10^5 


1.8 


350 


0, 


.71 


1.0 


5 


1 X 10^*' 


2.1 


350 


0, 


.71 


1.0 


6 


2 X lO^'' 


2.5 


350 


0, 


.71 


1.0 


7 


5 X lOif* 


2.8 


350 


0, 


.71 


1.0 


8 


1 X 10^'^ 


3.1 


350 


0, 


.71 


1.0 


9 


2 X 10^'' 


3.5 


350 


0, 


.71 


1.0 


10 


5 X lO^'' 


3.8 


350 


0, 


.71 


1.0 


11 


2.5 X 10^5 


1.5 


3500 


0, 


.71 


1.0 


12 


2.5 X 10^'^ 


2.5 


3500 


0, 


.71 


1.0 


13 


2.5 X IQi^ 


3.5 


3500 


0, 


.71 


1.0 


14 


5 X lO^** 


2.8 


350 


0, 


.95 


1.0 


15 


5 X 10^^ 


2.8 


350 


0, 


.50 


1.0 


16 


5 X 10^'' 


2.8 


350 


0, 


.05 


1.0 


17 


2.5 X IQis 


2.5 


350 


0, 


.71 


0.2 


18 


2.5 X lO^*^ 
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3.1. Temperature profiles 

The Figure [T] shows the temperature profile as a function of the Thomson optical depth resulting 
from constant density models for ten different ionization parameters (models 1-10 in Tabled]). In 
the Figure, the lower left curve corresponds to the least ionized case, and each consecutive to 
a higher value of F^. The values of log ^ are included next to the respective curves. Because 
the r grid is fixed with logarithmic spacing, the surface of the disk is in practice chosen to be 
Ttop = 10~^. The optical depth is measured from the surface towards the interior of the disk. Since 
the density is assumed to be constant, all these calculations are carried out along a distance of 
Az = Ar/axne ^ 1.5 x 10^° cm. 

These curves show very similar general behavior; the temperature is higher at the surface (due 
to the heating by the incident radiation) and decreases towards the interior of the disk, reaching 
a minimum temperature. For the high illumination cases (large values of £, and /or very close to 



the surface), the dominant mechanism is Compton heating and cooling (iKrolik et al.lll98ll ). with a 
rate given by: 



EF{E)dE - 4kT j F{E)dE 



(33) 



in the non relativistic limit. When this process dominates, the temperature approaches an asymp- 
totic value, the Compton temperature, given by the balance of the two terms in the previous 
equation: 

^ (34) 

where 

is the mean photon energy. From Equations (j34p and (|35p it is clear that Tq only depends on the 
shape of the spectrum. In particular, for the incident power law used in this paper (F = 2), the 
Compton temperature of the radiation is Tq = 2.8 x 10^ °K. Note that the temperature at surface 
for the case with log ^ = 3.8 (model 10) is indeed very close to Tq. 

In most cases, the gas remains at the high temperature limit throughout the region where the 
ionizing radiation field is unchanged from the incident field, thus forming a hot skin. The temper- 
ature gradually decreases in this hot skin as the bremstrahlung cooling becomes more important. 
As photons are removed from the Lyman continuum energies either by scattering or absorption 
and re-emission, at some point hydrogen recombines and the opacity grows rapidly, causing a drop 
in the temperature. Eventually the opacity becomes very large so that few photons are left in 
the 1-1000 Ry energy range, at which point the radiation field thermalizes to a nearly constant 
temperature (T ~ lO'' °K). 

It is worthwhile to mention that for the models with the two lowest ionization parameters 
(models 1 and 2), the opacity of the gas is large compared to the ionizing fiux even at the surface of 
the disk and therefore the gas reaches the low temperature limit very rapidly. Conversely, for the two 
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Fig. 1. — Temperature profiles for different illumination fluxes using constant density models (n = 
10^^ cm~^). The value of log ^ is shown next to each corresponding curve (models 1-10). 
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models with log ^ = 2.5 — 3.1 (models 6-8), the illumination is high enough to allow the photons to 
penetrate much deeper into the atmosphere, and the gas remains at high temperature (T ~ 10^ °K) 
even for large optical depths (r > 1). In fact, by looking the intermediate cases (models 3-8), one 
can clearly identify two temperature zones: the hot skin correspond to T > 10^ °K, while the 
cold region occurs for T < 10^ °K. Note also that although the transition between the hot and 
cold regions can be sudden, there is no thermal i nstability in these c alculations, since in all these 
cases the constant density restriction is applied (IKrolik et al.l Il98ll ) . Some of the models show 



small temperature inversions, such as the one in the curve correspoding to log ^ = 3.5 (model 9), 
that can be seen around tt ~ 6. This inversion occurs in the place within the slab where the 
temperature changes rapidly but also in a zone where the resolution in the optical depth is reduced 
due to la rger the step size in t he grid, which could lead to small scale thermal instabilities in the 
solution (iBuff &: McCravlll974l ). Despite this, there are no large scale thermal i nstabilities in thes e 
calculations, since in all these cases the constant density restriction is applied dXrolik et al.lll98lh . 



3.2. Reflected spectra 

Figure [2] shows the reflected spectra for each of the models shown in the previous section 
(models 1-10 in Table[l]), in the entire energy range covered in the calculations (1 eV to 210 keV). 
The spectra emerging from the top of the slab are plotted as solid curves, while the intrinsic disk 
flux is shown in the dotted curve. The dashed curve represents the X-rays incident at the top of 
the disk, assumed to be in the form of a power law with F = 2 and an exponential high energy 
cutoff at 200 keV (to improve clarity, only the incident power law for the lowest ionization case is 
shown). The corresponding values of log ^ are indicated above each curve, starting with the smaller 
value at the bottom, and increasing the ionization to the top. The curves are shifted to improve 
clarity. From bottom to top, each consecutive curve is rescaled by a factor increased by one order 
of magnitude with respect to the previous one (i.e., 1, 10, 100,..., 10^). Strong absorption profiles 
and edges are clearly seen in the neutral cases for log ^ =0.8, 1.1 and 1.5 (models 1-3), although 
it is important to mention that some small numerical problems occur in the thermal energy range 
of the spectrum for the lower flux case. For a few energy bins the reflected outgoing flux becomes 
negative, due to the large opacities for such energies (especially around 100 eV) and numerical 
errors in calculating the mean intensity. The redistribution of the photons due to Comptonization 
is evident above 10 keV (Compton bump), and in the smearing of the line profiles. For the models 
with the lowest illumination (models 1-3), there is a significant modification of the original power 
law continuum due to the large values of the photoelectric opacity for energies between 100 eV 
and 10 keV, where most of the strong absorption occurs. Nevertheless, none of these models shows 
a spectrum dominated by absorption. Even in the lowest ionization case (log ^ = 0.8), there are 
strong emission lines present through the whole energy range. This combination of emission and 
absorption features in the reflected spectra is a direct consequence of our accurate treatment of 
radiation transfer and the temperature gradient along the vertical direction of the disk discussed 
in the previous section. 
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The iron K line shows the effects of Compton scattering in models with log ^ > 1.5, since 
those are the cases for which a hot skin (T ~ 10^ °K) is present for at least a fraction of the total 
depth of the disk. By comparing the most ionized models, in particular those for log ^ = 2.8 and 
log ^ = 3.1 (models 7 and 8), one can see a very drastic change from a highly ionized to an almost 
featureless spectrum, and even fewer features are seen in the two models for log ^ = 3.5 and 3.8 
(models 9 and 10). This is due to the fact that in these cases the gas is always at high temperature 
within the range of our calculations. Despite this, emission from highly ionized iron K lines is still 
apparent in the reflected spectrum. 

Figure [3] shows the same results as Figure [21 but in the 2-10 keV energy range. Here it is 
possible to see the structure of the iron inner shell transitions in detail. The Compton shoulder 
is apparent in the first three curves from bottom to top (least ionized cases). Changes in the Fe 
Ka emission line are distinguishable as the ionization parameter increases. The neutral K/? line 
is only evident in the three least ionized cases (models 1-3), as well as the Fe absorption edge at 
~ 7.5 keV. This figure also shows a more complex structure of the iron line in intermediate cases of 
ionization: line positions and intensities vary between mostly neutral emission at 6.4 and 7.2 keV 
for Ka and K/3, to a very ionized profile with lines around 6.5-7.0 keV and 7.8-8.2 keV for the same 
transitions. 

We have calculated the equivalent widths for the Fe Ka lines in all these spectra. Because of the 
large deviations of the reprocessed continuum from the incident power law in some of these models 
(particularly for low ^), in order to calculate the equivalent width we define a local continuum by 
interpolating a straight line between two points in the region around the line, specifically between 
6 and 7 keV. The integration is performed within this range as well, such that only Ka emission is 
taking into account. 

The resulting equivalent widths for the Fe Ka emission line vary between ~ 400 — 800 eV for 
the cases with 1 < log.^ < 3 (left panel in Figure [3]). For higher values of the ionization parameter 
the equivalent widths decrease very rapidly, to approximately 40 eV for the most ionized case 
(log ^ = 3.8, model 10). Since the gas is more ionized for high illumination, the emission of the line 
becomes exclusively due to H- and He-like iron ions, which combined with the extreme Compton 



scattering of the photons results in the reduction of the ec 



resembles the X-ray Baldwin effect (llwasawa &: Taniguchi 



active galactic nuclei spectra (e.g., iPage et al.ll2004l ). 



uivalent width of the line. This tendency 



19931 ). which have been observed in many 



3.3. Spectral features 

In order to provide a finer analysis of the refiected spectra, we perform calculations increasing 
the energy resolution by an order of magnitude, i.e., using 5 x 10^ energy grid-points {TZ 3500). 
Since these calculations are more expensive in terms of CPU-time, we only show three representative 
cases. These calculations are shown in Figures SHU and correspond to models 11, 12 and 13 
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Fig. 2. — Solid curves: Reflected spectra for tlie models presented in Figured) The value of log ^ is 
shown next to each corresponding curve (models 1-10). The curves are shifted by arbitrary factors 
for clarity. These are, from bottom to top: 1, 10, 100,. ..,10^. Dashed curve: incident ionizing 
spectrum (power law), corresponding to the lowest ionization case (log ^ = 0.8). Dotted curve: 
intrinsic disk flux (black body), common to all the models. 
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Fig. 3. — Reflected spectra as shown in Figure [2] in the 2-10 keV region. The value of log ^ is 
shown next to each corresponding curve (models 1-10). The curves are shifted by arbitrary factors 
for clarity. These are, from bottom to top in each panel: 1, 10, 100, 10^, 10^. 
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in Table [TJ These spectra (Figures SHI]), are plotted in their physical units and without any 
renormalization or rescaling, and besides the much higher energy resolution, all the other input 
parameters in the models are the same as the ones used for the models present previously in 
Figures [1] and [2j 

Figure m shows the reflected spectra corresponding to log ^ = 1.5 (model 11), in the 0.5-10 keV 
energy range, while using a resolving power of 7^ 3500. The most prominent and relevant emission 
lines are identified in the Figure. In general, the spectrum is dominated by many absorption edges 
plus radiative recombination continua (RRC). RRCs occurs when a free electron is captured by an 
ion into an unoccupied orbit; if the electron carries more energy than it needs to be bound to the 
ion, the excess will be radiated as a photon. Emission K lines from H- and He-like nitrogen, oxygen, 
neon and magnesium are present in the region around and below 2 keV. At higher energies, there is 
clear emission from Si xii-xxiv and S xv-xvi. The iron Ka and K/3 components are clearly visible 
at ~ 6.4 keV and ~ 7.1 keV, respectively, as expected from mostly neutral emission. Many lines 
are blend together in a small region (about 50 eV for each component), but the emission is mainly 
due to Fe xiii up to Fe xvi. The Compton shoulder in both components is also evident, as well 
as a very marked absorption edge at ~ 7.5 keV blended with a complex structure of absorption 
profiles. These absorption features correspond to resonances given by transitions from the ground 
state to the Is — np autoionizing states, where n = 3, ...,30; for our particular atomic data set. 
These are K-vacancy states that can decay either radiat ively or via Auger spectactor channels. 



being the latter dominant for n > 3 (jPalmeri et al.ll2002l ). The resonance structure seen in this 



spectrum covers the 7.1-8. keV energ y range, thus i t is ma inly due to second row Fe ions, which 



resemble the opacity curves shown by iPalmeri et al.l (l2002l ). as well as the photionized spectra in 



Kallman et al.l (j2004l ). for similar values of the ionization parameter considered here. The energy 



position of the absorptioi i edge indicates that is likely to be produced by Fe xili-xiv, according to 



the energies presented in lKallman et al.l ( 2004 ) (see in particular their Figure 2 



Figure [5] shows the resulting spectra when log ^ = 2.5 (model 12), again for a high resolution 
energy grid {TZ ~ 3500). There is a significant change in the overall continuum, particularly in the 
1-10 keV region, in comparison with the previous case (log ^ = 1.5). This is because in the lower 
ionization case the photoelectric opacity dominates in this energy band and changes the original 
incident power law shape, while the continuum spectrum in Figure [5] still shows the original slope 
(r = 2). Lower opacity also allows the efficient propagation of lines (since photons are able to 
escape), and the emission from heavier elements (Ar xvii, Ca xix and Ca xx), since the gas is 
more ionized. In general, this spectrum shows more emission lines and weaker absorption edges. 
However, some RRCs can still be identified, especially around the sulfur and calcium lines. The K 
lines from nitrogen, oxygen and neon ions are again the strongest emission features in the low energy 
part of the spectrum. The high energy region changes noticeably with respect to the previous case 
with lower ionization. The strongest Fe Ka emission in this case is due to Fe xxv ions at ~ 6.7 
keV. However, there is a rich structure of emission lines with centroid energies that cover from 
~ 6.4 keV to 6.7 keV produced by ions at lower ionization stages (Fe xv-xviii). The Fe xxvi 
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Fig. 4. — Reflected spectra for log ^ = 1.5 while using a resolving power of 7^ ~ 3500 (model 
11). All the other input parameters are the same used before (as in Figures [1] and [2]) , although no 
rescaling is applied. The strongest emission lines are labeled. 
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Kq emission line is clearly observable at 6.965 keV. All these lines are superimposed on a very 
broad, Comptonized profile. The K/3 component is replaced by a few distinguishable emission lines 
mainly due to H- and He-like iron. As in the previous case, a rich structure of absorption profiles 
is noticeable between 7.2 and 8.2 keV. A small absorption edge can be seen just below the RRC 
located at ~ 9 keV. 

It is important to notice that the O viii Ka emission line shows the same kind of broadening 
as the Fe K line. In fact, the other oxygen lines also show broadening in some degree, as well 
as those coming from neon. This is somewhat unexpected, since photons at those energies would 
need many scatterings in order to produce such broadening if the gas were cold. However, since 
the temperature is high (T ~ 10^ °K) even at r ~ 1, broadening also takes place through the 
dependence of the Gaussian kernel on the temperature (see § 12. ip . This also means that the H-like 
oxygen lines are being produced efficiently over a large range of optical depths, which can be verified 
by looking at higher ionization parameters. Figure [6] shows the last model with a resolving power 
of 7^ 3500, resulting from a gas at log ^ = 3.5 (model 13). Almost no absorption features can 
be seen in this spectrum, and the emission comes mostly from recombination of fully ionized ions. 
Both iron and oxygen Ka line profiles are very broad and Comptonized. There are no apparent 
signs of broad components in the weaker lines, such as those from neon, although this is likely 
due to the fact that the redistribution is so extreme that it completely smears the profile over the 
continuum. Therefore, only those photons produced very close to the surface come out of the slab 
unscattered. The emission in the iron K region is exclusively due to H- and He- like ions. There is 
strong emission from the Rydberg series of Fe xxv, namely transitions from Is 4p, 5p and 6p to 
the ground state Is^, at 8.62, 8.83 and 8.94 keV, respectively. 



3.4. Anisotropy: incident and view^ing angles 

The boundary condition expressed in Equation (I14p depends, in general, on the angle with 
respect to the normal on which the radiation is incident. Isotropic illumination implies that the 
source of radiation is extended and very close to the region where the calculation takes place. 
On the other hand, locating the source of the external X-rays at a specific point far away of the 
accretion disk constrains the angle at which radiation will penetrate in the atmosphere. To account 
for such a situation, the right hand side of Equation (I14p can be expressed as 

Iinc = lo^ifJ.- fio), (36) 

where /io is the cosine of the incidence angle. The lower panel of Figure [7] shows the refiected spectra 
for a constant density model in which log ^ = 2.8 for three different incidence angles (models 14-16 
in Table [1]). The red curve corresponds to = 0.95 (normal incidence), the green curve correspond 
to ;Uo = 0.5, and the blue curve to /io = 0.05 (grazing incidence). The upper panel of the Figure 
shows the corresponding temperature profiles as a function of the Thomson optical depth, obtained 
for each case. At the surface of the disk (r = 0), the temperature increases with the incidence 
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Fig. 6. — Reflected spectra for log ^ = 3.5 while using a resolving power of 7^ ~ 3500 (model 
13). All the other input parameters are the same used before (as in Figures [1] and [2]) , although no 
rescaling is applied. The strongest emission lines are labeled. 
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Fig. 7. — Temperature profiles (upper panel) and reflected spectra (lower panel), resulting from a 
constant density model with log ^ = 2.8 for three different incidence angles (models 14-16). In the 
two panels, the red curves corresponds to = 0.95 (0 ~ 0°, normal incidence); the green curves to 
^0 = 0.5 (or 6 = 60°); and the blue to /^o = 0.05 {0 90°, grazing incidence). 
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Fig. 8. — Outgoing intensities in the 2-10 keV energy range for constant density models from 3 
different ionization parameters (models 2, 5 and 8), as viewed at three different angles with respect 
to the normal. In all cases, the red curves corresponds to // = 0.95, the green curves to = 0.5, 
and the blue curves to /i = 0.05. 
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angle, such that the grazing incidence shows the highest temperature of the three cases. However, 
when the illumination occurs at normal incidence, the radiation ionizes deeper regions in the slab 
and there is a larger amount of hot material for those cases. In this model, the physical quantity 
that describes the amount of illumination is the net X-ray flux Fx, in units of energy per unit area 
per unit of time. Since the flux is defined as the second moment of the intensity, (Equation I24p . 
therefore Fx = ^ J lincfJ-dfJ-, and using Equation (f36]l 

line = — - fio), (37) 

which means that for the same value of (or the same ionization parameter varying the 
incidence angle fiQ effectively varies the intensity of the radiation incident at the surface. This will, 
of course, produce more heating and raise the temperature and ionization of the top layers in the 
disk. However, the angular dependence of the radiation field ensures that in the cases of normal 
incidence the illuminating radiation reaches deeper regions in the atmosphere than those of grazing 
angles, as is expected. Although the effects on the reflected spectra are not obvious, the spectrum 
for the grazing incidence (blue curve) is stronger than the other two, mimicking a case of a higher 
illumination (note that none of the curves have been offset). Nevertheless, the emission lines are, 
in general, very similar in all cases. It is also interesting to see how the temperature profiles of 
the three models converge at large optical depths, since the radiation fields become more isotropic 
after many scatterings. 

The anisotropy of the reflected radiation field can also be investigated by looking its angular 
distribution through different viewing angles. Figure [8] contains the reflected spectra for 3 different 
ionization parameters (indicated next to each case), as is observed from 3 different angles with 
respect to the normal, in the 2-10 keV energy range. For consistency the colors represent the same 
angles as in Figure[Tl i.e., the red curves indicate fi = 0.95 (face-on), the green curves n = 0.5, while 
the blue curves correspond to ^ = 0.05 (edge-on). Note that no rescaling or normalization is applied 
to these curves, therefore the differences are exclusively due to the differences in the ionizing fluxes 
and on the viewing angles. Moreover, the physical quantity plotted here is the outgoing specific 
intensity /(O, E) instead of the flux (as the other Figures), which is an angle averaged quantity. 
In fact, these 3 cases are those presented in Figures [2] and [3] (models 2, 5 and 8 in Tabled]), and 
thus their temperature proflles correspond to those shown in Figured] In general, the reprocessed 
features are stronger when th e disk is observed face-on than when the viewing angle is parallel to 
the surface. Nayakshin et al. (j200ol ) found the same tendency in their hydrostatic models, which 
they explained to be a consequence of an effective Thomson depth that changes when the slab is 
observed at different angles, according to the geometrical projection Te// = t/h. However, their 
results suggest that the reflection features almost disappear for small values of /i, and that no 
visible iron could been detected in such cases. The results shown in Figure [8] do not completely 
agree with this, since the Fe K lines are strong even for the high ionization case. Nevertheless, one 
must take this comparison wit h care, since all the res ults presented in this paper correspond to 
constant density models, while iNavakshin et al.l (|2000l ) calculations were done under hydrostatic 
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equilibrium. In any case, the trend in the emerging spectra shown in Figur es [7] a nd El is similar to 



previo us models in the literature, such as those by iNavakshin et al.l (|200d ) and iBallantvne et al 
(|200lh . 



3.5. Iron abundance 

We have also studied the effects of the iron abundance on the ionization of the atmosphere 
and the reprocessed spectra. The upper panel of Figure shows the temperature profiles along the 
vertical direction for a gas at log ^ = 2.5 when the iron abundance is assumed to be 0.2, 1, 2, 5 
and 10 times its solar value (models 17, 6 a nd 18-20, respec t ively) . The solar atomic abundances 



used in all our calculations are those from iGrevesse et al.l (|l996l ). and in particular for iron is 
3.16 X 10^^ (with respect to the hydrogen). The curves towards the left of the plot correspond to 
higher values of the iron abundance, which means that the gas is effectively cooler than the solar 
and sub-solar cases. This shows that the iron produces a net cooling as its abundance is increased. 
The transition between the hot and cold regions in the gas occurs at lower optical depths for the 
super-solar abundance models, decreasing the thickness of the hot skin. However, all the models 
converge to the same temperature in the cold region of the disk (T 2.5 x 10^ °K). 

The corresponding reflected spectra are shown in the lower panel of Figure [H These curves are 
shifted by arbitrary factors for clarity, which are, from bottom to top: 10"^, 10°, 10^, 10^, and 10^. 
The iron emission is enhanced when the abundance is increased, as can be clearly seen at ~ 6 keV, 
~ 1 keV and ~ 0.1 keV for the K-, L- and M-shell transitions, respectively. It is also evident that 
the continuum is highly modified in the 1-40 keV energy region, which is where the larger opacity 
takes place. However, there are no significant modifications at higher energies (i.e., if these spectra 
are placed on the top of each other, there are no differences in the region above 40 keV). 

The suppression of the continuum combined with the enhancement of the K-shell emission 
lines as the iron abundance is increased shows a significant impact on the line equivalent width. 
When solar abundance is assumed, we found an equivalent width for the Fe Ka of ^ 700 eV. This 
value is increased to about 1.1, 2.3 and 3.9 keV when iron is chosen to be twice, five and ten times 
more abundant, respectively. Conversely, this value decreases to ~ 203 eV for the model when 
Ape = 0.2. 



3.6. Comparison with previous models 



Taking into account that our main contribution to the already existing models of accretion 
disk is the inclusion of the most complete atomic data available (particularly important for iron), 
we compare our results for constant density accretion disks models with those included in reflion 
(jRoss &: FabianI 120051 ) . Such a comparison is shown in Figure [TOl where the reflected spectra for 
three different ionization parameters are shown (log ^ = 1.5,2.5 and 3.5). In the figure, the black 
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Fig. 9. — Upper panel: temperature profiles from a constant density model with log ^ = 2.5 and 
different iron abundances with respect to the solar values. From right to left, each curve correspond 
to: AFe=0.2, 1, 2, 5 and 10 (models 17, 6 and 18-20, respectively). Lower panel: reflected spectra 
for the same models. The curves are shifted by arbitrary factors for clarity. These are, from bottom 
to top: 10~^, 10°, 10^, 10^, and 10^. The values of the iron abundance are shown at the top of each 
curve. 
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Fig. 10. — Reflected spectra from constant density models for three different ionization parameters 
(models 11-13), plotte d in the 2-10 keV energy range. The red curves are the models from reflion 
(|E.oss fc Fabianll2n0,'^ ). for the same input parameters. 
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curves correspond to our calculations (models 11-13), while the red curves are the models from 
REFLION with the same input parameters. It is clear from this comparison that these two models 
show important differences. The Fe K line profile around 6.4 keV from our model is shown to be 
much more intense but with similar widths, and the absorption edge at ~ 7.5 keV to be quite strong 
(especially for the lower ionization). In fact, the equivalent widths of the iron Ka emission line are 
similar in both models for log = 1.5, our models predict a value of 610 eV while reflion predicts 
793 eV. However, there are larger discrepancies in the higher ionization cases: our models predict 
equivalent widths of 695 eV and 92 eV for log = 2.5 and 3.5, respectively; while reflion models 
show values of 473 eV and 49 eV for the same cases, respectively. As before, these equivalent 
widths were calculated using a local continuum which is defined by interpolating a straight line 
between two points. In this case the integration was performed between 6 and 7 keV in order to 
take into account only the contribution from the Ka emission (since any Fe K/3 emission will occur 
at energies above 7 keV). Furthermore, the same integration region is applied to all the models in 
order to provide a consistent quantity to compare models with different ionization parameters. 

However, we have also calculated the equivalent widths when the integration region is modified 
for each model such that the entire broadened line is included. The new integration regions are: 
6.1-6.6 keV, 5.9-7.3 keV, and 5.2-7.8 keV for models with log ^ = 1.5,2.5 and 3.5, respectively. 
The resulting equivalent withs from our models are 569 eV, 829 eV and 255 eV; while reflion 
predicts 700 eV, 813 eV and 420 eV for the same models, respectively. Clearly, there is a better 
agreement when the entire line profile is taken into account, especially for the two lowest ionization 
cases. Although both set of models share the same input parameters, in reflion the illuminating 
spectrum extends up to 1 MeV with an exponential cutoff placed at 300 keV, while our models 
cover up to 210 keV, with the exponential cutoff placed at 200 keV. This could be responsible 
for some of the differences seen here, in particular by altering the broad component of the Fe K 
emission which is highly affected by the Comptonization of the photons in the gas. In general, our 
models predicts stronger Fe Ka emission, a more complex structure for the Fe absorption edge, 
and the presence of the Fe K/3 emission line at 7.2 keV, not seen in the reflion models. 



Additionally, we have also compared our models with those presented by iNayakshin et al 



(|2000l ). In many aspects these models are the most similar to the ones presented in this paper, 
in particular since they have also used the code xstar to solve the ionization structure of the 
gas. By means of visual comparison of their figures we found a good general agreement in the 



five constant density models included in iNavakshin et al.l (j200d ) . However, our models predicts 



higher temperatures in the regions close to the surface of the disk (at small optical depths), and 
lower temperatures at large optical depths; which affects the ionization of the gas and the final 
reprocessed spectrum. This is likely due to the fact that their models only explicitly treat a smaller 
optical depth (rxhomson = 4) than ours, and the radiation transmitted through this region is treated 
approximately. 
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4. Conclusions 

Answers to outstanding questions concerning accretion disks will not come alone, but with the 
need to reproduce and predict new observations and to identify spectral profiles at high resolution. 
In this paper we have presented new models for the structure of X-ray illuminated accretion disks 
and their reflected spectra, assuming constant density along the vertical direction. These models 
include the most recent and complete atomic data for the iron isonuclear sequence. The energy 
resolution used in the reflected spectra exceeds other models previously published as well as the 
resolution of the detectors on current X-ray observatories ( Chandra, XMM Newton, Suzaku) , and it 
is comparable to the expected resolving power of the forthcoming International X-Ray Observatory 
(IXO). 

In models with intermediate values of the ionization parameter (1.5 < log^ < 3.1), the structure 
of the gas often displays a two temperature regime: a hot skin (T > 10^ °K) close to the surface 
where the Compton heating and cooling dominates, and a cold region (T < 10^ °K) where the 
photoelectric opacity quickly thermalizes the radiation fields. The thickness of the hot skin increases 
with the illumination, and although these solutions are thermally stable the transition between the 
hot and cold regions can be sudden. The emerging spectra corresponding to these models show 
a clear combination of absorption and emission features, particularly for log ^ < 2.5. For higher 
illumination cases the effects of Compton scattering become more evident, even in the low-energy 
part of the spectrum {E < 100 eV). In the high energy part of the spectrum Compton scattering 
partially smears the profile of the iron K line. The Compton hump above 10 keV is clearly visible 
in all the models considered here. 

The equivalent width of the Fe Ka emission line varies between ~ 400-800 eV for models with 
1 < log^ < 3. Once the ionization parameter is larger than 10^, the equivalent width decreases 
drastically. The lowest value of the equivalent width is ~ 20 eV, corresponding to the highest 
illumination case (log ^ = 3.8). This behavior resembles the X-ray Baldwin effect, recently observed 
in several AGNs. A comparison with other models previously published such as REFLION shown 
important differences in the structure of the iron K lines. In particular, the equivalent widths in 
those models can differ by a factor of 2 when compared to the ones shown in this paper. 

In general, our simulations show that the K-shell atomic data is crucial to properly model the 

structure and profile of the iron lines. These results also suggest that the line emission from Fe 
ions in different ionization stages and Comptonization of high energy photons by cold electrons can 
be responsible for significant line broadening. These processes need to be taken into account since 
they can be mistaken for relativistic effects, especially in cases when the gas is partially or high 
ionized. 

Because the state of the gas in the accretion disk depends on the radiation field at each point, 
and the ionization balance required for a realistic definition of the source function in the radiation 
transfer depends on both the temperature and density of the gas, a self-consistent approach must 
be taken into account in order to solve the hydrostatic equilibrium equation instead of assuming 
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constant density. This will be considered as the next step of this project. Future work will also 
consid er the inclusion of new K-s hell atomic data recentl y calculated for the Ne, Mg , Si, S, Ar and 



Ca by iPalmeri et al.l ( 2008all and IWitthoeft et al.l (|2009l ) , as well as for the nickel (jPalmeri et al 



2008bl ). and nitrogen (jGarcia et al.ll2009l ) isonuclear sequences. 
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